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The rotation of the bodies and the eccentricity of the orbit have significant effects on the emitted gravita- 
tional radiation of binary systems. This work focuses on the evaluation of the gravitational wave polarization 
■ states for spinning compact binaries. We consider binaries on eccentric orbits and the spin-orbit interaction 

up to the 1.5 post-Newtonian order in a way which is independent of the parameterization of the orbit. The 
equations of motion for angular variables are included. The formal expressions of the polarization states are 
given with the inclusion of higher order corrections to the waveform. 

1 Introduction 
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With the construction of the ground based gravitational wave observatories, i.e. LIGO [1 , VIRGO 2 , GEO600 
[3] and TAMA [4], the investigation of gravitational radiation has been receiving new impetus in the last two 
decades. The description of the dynamics of different astrophysical sources and the detectable signal of the 
emitted gravitational waves has become one of the main lines of research in general relativity and astrophysics. 

Among the plausible sources of gravitational waves which can be detected by these observatories are binary 
systems of compact objects. Having the appropriate amplitude and frequency properties for their detection, 
fNj ■ ground based detectors are sensitive mainly for the radiation emitted in the merger era. These kinds of signals 
are described both by numerical simulations 5 and analytic approaches 6,7. Beyond the Lazarus approach 
\Q • [SI [3] which combines numerical techniques and perturbation theory to extract approximate waveforms, recent 
progress in numerical relativity have made it possible to simulate the last stages of inspiral, merger and ringdown 
in full general relativity [TU1CEI]. The comparison of numerical and analytical models shows good agreement in 
the gravitational wave phase evolution in the two approaches [12l [13] and helps us to check the accuracy of the 
analytical models and to discriminate between them [14] . 

Furthermore, for the detection of gravitational waves the description of the inspiral phase of the coalescence 
is also relevant. The main purpose of the determination of the detectable gravitational wave signal in the 
inspiral era is the LISA project [131 [H], which aims to install a detector system into space. Because of their 
different size and technology LISA and the other space based detectors will be most sensitive in the frequency 
range of the emitted radiation of inspiralling compact binaries. Because of the short duration of merger it is 
important to determine the direction of the source as soon as possible [T7]. An additional advantage of this 
description is to give analytic initial conditions for the simulations modeling the coalescence. 

During inspiral the dynamics of the source and the emitted gravitational radiation are described by two 
highly related methods. The post-Newtonian (PN) approximation [T51 US] is used to investigate the dynamics 
and the motion of the source. The essence of this approach is the neglection of ultrarelativistic and extreme 
gravitational effects with the introduction of a series expansion, where the expansion parameter is e ~ v 2 ~ M/r 
with the typical velocity v, size r and total mass M of the system. This approximation can be used until the 
binary reaches the innermost stable circular orbit [201 HI] • To describe the emitted gravitational waves the post- 
Minkowskian approximation is applied. In the outer region of the source and the near zone of the radiation the 
two expansions can be matched together [22] . 

The dynamics of the binary is well described in the post-Newtonian approximation. The conservative part 
of the dynamics is known to high post-Newtonian orders [23 . The equations of motion with the inclusion of 
spin-orbit effects and the spin precession equations are given in |24j beyond the leading-order contributions, up 
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to 2.5 PN order. The radial motion and the parameterization of the orbit are known up to 3PN order [25]. 
However, these results are valid under the assumptions that the rotation of one or both of the objects and 
radiation reaction are negligible. For spinning compact binaries the dynamics and radiation reaction effects are 
discussed up to 1.5 PN order [26]. The full description of motion which contains the evolution of the angular 
variables is known only in special cases, namely, in the test particle limit [27] , for binaries with one spinning 
components [28] up to 1.5 PN order, or with the neglection of spin effects up to 3PN order [2"3] , 

As a result of the gravitational wave research the transverse-traceless part of the emitted radiation field is 
given formally in terms of the main dynamical quantities of the binary system (the separation, relative velocity 
and spin vectors), the mass of the bodies and the direction of the line of sight. The results including spin 
effects are valid up to 1.5-2 PN order [30, 31 . Higher order contributions to the waveform are known in the 
non-spinning and circular orbit case up to 3-3.5 PN order [23]. A detailed investigation of the relativistic 
post-Newtonian corrections with the inclusion of higher multipole moments for eccentric orbits is given in [32] . 

In our previous works we have given the description of the dynamics and a method to evaluate the detectable 
waveform of the emitted gravitational radiation for binaries with one spinning component and in the test 
particle limit [28] [33] . Since the detectable gravitational wave signal can be unambiguously decomposed into 
the independent polarization states h + and h x as 

h(t) = F+h+(t) + F x h x (t) , (1) 

where F+ and F x are the beam pattern functions [34] , the main goal of our work is to construct a consistent 
method to determine the explicit form of these polarization states. 

In the present paper we generalize our former results [28j to discuss binary systems with components of 
arbitrary masses nij and ni2 and spins Si and S2. Our method and formulae are valid up to 1.5 PN order and 
applicable to eccentric orbits independently from the parameterization of the orbit. 

This paper is organized as follows. In Sec. II we present the description of the elements of the dynamics 
used in the projection method. In Subsection II A we introduce the main coordinate systems and the angular 
variables. Subsection II B contains the solution of the spin precession equations. In Subsection II C we discuss 
the equations for the components of the relative velocity vector and the equations of motion for the angular 
variables. Sec. Ill contains the main results of our paper, namely, the projection method and the expressions of 
the polarization states of the detectable gravitational waves. First we give the components of the vector triad 
used in the projection in Subsection III A. In Subsection III B we give the formal expressions of the polarization 
states in terms of the main vector components and constants of motion. Here we outline the main steps of the 
method which leads to the polarization states. In Sec. IV we discuss the circular orbit limit, where the explicit 
time dependence of the dynamical variables and the waveform can be given. In the last section we summarize 
our results and give a detailed outlook. 

2 The motion of the binary system 

To describe the dynamics of a binary system and the polarization states of the emitted gravitational waves one 
needs two different coordinate systems. The equations for the orbital elements can be easily evaluated in the 
invariant system which does not change in time, moreover, the projection method leads to simple expressions 
in the comoving system. 

The following description of the evolution of the binary provides us the equations of motion for the angular 
variables describing the direction of the separation vector r and gives short expressions for the different com- 
ponents of the relative velocity vector v. All these quantities depend on two scalar and one vector constants of 
the motion, i.e. the energy E of the system, the length of the angular momentum vector L = Ltv + Iipn + L50, 
where Ljy = A*r x v is the Newtonian angular momentum, fi — m,\milM is the reduced mass, M = mi + 
is the total mass of the binary, and Lpjv and ~Lso are the 1 PN and spin-orbit contributions, respectively. The 
constant vector of motion is the total angular momentum J = L + Si + S2. All these quantities and the length 
of the spin vectors are constant up to 2 PN order. 

2.1 Angular variables 

The dynamics of the binary and the evolution of the orbital elements are most conveniently described in the 
invariant coordinate system. The z axis of this system is fixed to the direction of the total angular momentum 
vector J. There is another constant vector in our description, the direction of the line of sight N. To simplify 
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the calculations the x and y axes of the invariant system are chosen in a way that N = (sin 7, 0, cos 7) in this 
system with the constant angle 7 of J and N. 

The projection method results in simple expressions for the polarization states if we express the components 
of the main vector quantities and the transverse traceless part hj, T of the radiation field in the comoving 
coordinate system. The x and z axes of this coordinate system are aligned with the separation vector r and the 
direction of the Newtonian angular momentum L^v, respectively. 

The transformation between these coordinate systems is characterized by Euler-angles [35] . The vector u in 
the comoving system is parameterized as 

u' = fl,($)fl l (t)J?,($)u (2) 

in the invariant system, where 1 is the angle between J and Ljy, $ describes the precession of the orbit and \P 
is the polar angle on the orbital plane. Using this transformation the separation vector in the invariant system 
becomes 

cos $ cos ^ — cos 1 sin $ sin iff 
r I sin <!> cos \& + cos 1 cos $ sin \I/ | (3) 
sin 1 sin if> 

In the comoving system the components of the relative velocity vector are expressed by Euler-angles, 

( ^ \ 

v = r (cos + \IM . (4) 

V / 

For further convenience we decompose the relative velocity vector as 

v = U||n + tiiDi , (5) 

where n = r/r is the direction of the separation vector and m is the unit vector in the direction of the y axis 
of the comoving system. Moreover, to separate terms which contribute to different PN orders and the rotation 
of the bodies we decompose the main quantities as follows 

r = r N + r PN + r s , v\\ = v\\ N + v\\ PN + v\\ s , v± = v ±N + v ±PN + v±s , 

* = *AT + fpAT + *S , $ = $jy + *PJV + $S , ' = ijV + l-PN + , (6) 

where r denotes the length of the separation vector. 

During the calculation of the polarization states we do not need the expressions for all the Euler angles 
separately. As it will be shown below only the combinations T = $ + <f, 15 sin at and tscos$jv appear. To 
determine the equations for these quantities we note that ln = 0, see the Appendix of Since the 1 PN 

corrections do not change the direction of the orbital plane l P n — 0, and hence 1 = l$- 

2.2 Spin-precession 

The first step in describing the dynamics of the binary system is the examination of the spin-precession equations, 

Si = ^(^L w -S 2 + i(r.S 2 )r)x Sl , (7) 
1 ^4 + 3C2 T 3 



S 2 = - ( -V^Ljv-Si + ^r-SOr) x S 2 , (8) 



where £i — C2 1 — m^frnx. 

Since we are interested in the lowest order contributions (these contributions give rise to 1.5 PN effects in 
the description of motion) the spin-spin interaction which appears at higher order is neglected and hence 

S.^JxS, (9) 

for i = 1,2. 
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The components of the spin vectors are expressed with the magnitude Si and the spherical angles on and 
ft. The evolution of these angles is governed by 



- sin ft sin o,ft + cos ft cos otiCHi 
Si I cos Pi sin a,ft + sin Pi cos ciiOn 
sin a^di 



JS- 

= (4 + 3Q)—-l \ cosftsma-i 




(10) 



where the lengths of the spin vectors change only under radiation effects over 2PN order. In this way we get 
the following equations for the spherical angles 



0, ft = (4 + 36)^ 



(11) 



In the Newtonian case Tn governs the angular dynamics of the system. With the use of Eq. (fTTj) we relate 
the relative PN order of the angular variables to it, 



ft _ (4 + 3C l )J/(2r%) L/(r%) _ f i 
T N L/{fir 2 N ) " L/(prff) r N 



(12) 



Since we require that the integration of the differential equations does not change the relative PN order of the 
different quantities the same estimate is right for the angular variables, henceforth ft/Yjv ~ e- This statement 
indicates that we can introduce the decomposition ft = ftjv + PiPN, where the relative order of these angles 
are ftjv/T n ~ 1 and PiPN /^n ~ £■ Here ftjv are constant and 
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2.3 The equations of motion 

If one combines the results of [55] and [31] the length and the components of the relative velocity vector can be 
written as 
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where 77 = fi/M and we use the notations <r = C1S1 + C2S2 and S = Si + S2. Since v 2 _ = v 2 — r 2 we have 
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To derive the equation for the angle T we use the form of the relative velocity vector given in Eq. ^ and 
that l = bg. From the second component we have v± — rT, and hence 
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To evaluate the remaining equations the components of the total angular momentum J = Lat+Lpjv+Lso+S 
are determined, where 
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is the 1 PN correction of the angular momentum vector and 

Lso = V j -s-[r x (r x (2S + cr))] - I[ v x (v X cr)] 



(18) 



is the spin-orbit contribution. In this way the components of J are 

Ll s sin$jv + [Ai sin^ + Bi cos^ + (d + 1) cos (3iN + PiPN sin/3ijv] Si sinai 
t=i 

2 

-its cos$/v + ^ i B i sm d ~ Ai cos£ + (Ci + 1) s'mfiiN - (3 iPN cos f3 iN ] Si sin a, ; 

2 

L + ^(C 4 + l)S l cosa t 
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where we have introduced the notation = 2T/v — and 
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Since the total angular momentum vector is J = (0, 0, J) in the invariant system these equations give rise 
to the identities 

2 2 

Si cos PiN sin aj = , 5,; sin /3jjy sin a>i = (21) 



in 0.5 PN order, and 

ts sin <£> jv 



i5 cos $ 



A' 



1 2 

- — 2J sin & + Bi cos f j + d cos ftjv + ftpiv sin /3w) Si sin ctj 

i=l 

I 2 

- ^ {Bi sin & - ^ cos + Ci sin (3 iN - f3 iPN cos /3 iA r) Si sin a< 



(22) 



in 1.5 PN order. As it is not necessary to evaluate ts and $at separately with Eqs. (|22p we have all the 
expressions required for the description of motion of the binary system and the emitted gravitational radiation. 



3 The polarization states 
3.1 The projection 

The polarization states of the detectable gravitational waves are determined by the projections of the transverse- 
traceless tensor hl^ T representing metric perturbations, 



h+ = ^{PiPj - qiqj)hji T , h x = -(ptqj + qiPj)h% T , 



(23) 



where p is a unit vector in the orbital plane perpendicular to the direction of the line of sight N, and q = N x p. 

Since we have chosen the comoving system to describe this projection we determine the components of N, 
p and q in this system. With the use of the transformation law @ N has the following form in the comoving 
system 

sin 7 cos — [(Tpjv + Ys) sin 7 sin Tjv — cos 7 sin <&at<,s] 
N= ( — sin7sinTAr — [(Tp/v + T5) sin7CosTjv — cos7COS$Art5] J . (24) 

cos 7 + SU17 sin 



Moreover, the conditions for p determine its components as 

sin Tat + (T PN + T s - ctgj cos $/vis) cos Tjv 
p = ( cos Tjv - (Tp/v + T s - ctg 7 cos $nl s ) sin Tjv 





(25) 



From Eqs. (|24|) and (|25|) it follows that for the description of the polarization states the relevant angular 
variables are T and the combinations (.scos^/v and tssin$/v. 
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3.2 The expressions for the polarization states 

To separate the different contributions the transverse-traceless part of the radiation field h l ^ T can be decomposed 
in the post-Newtonian approximation as [30) 

2\x 
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(26) 



where D is the distance between the observer and the source and we have collected all the terms which are 
relevant up to 1.5 PN order. Q u denotes the quadrupole (or Newtonian) term, P°- 5 Q* J ', PQ W and P 1 - 5 ^ 1 - 7 are 
higher order relativistic corrections and PQ % so anc ^ P Q % so are ^ ne spin-orbit terms. The detailed expressions 
for these contributions are given in [3U1 [3T] • 

We choose a similar decomposition for the polarization states h + and h x , 
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It is worth to mention that this decomposition is formal in the sense that higher order terms in the description 
of motion can cause higher order contributions in h+ and h+ ' 5 . 

With the use of the decomposition (J5]) of the relative velocity vector we evaluate the projections in Eq. (|23|) 
and obtain the following formulae for the "plus" polarization state 
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where 8m = m 2 — m\ and A = M(S 2 /m 2 — Si/mi). 

The main steps of the method, which results in the detectable gravitational wave signal, are the following. 
With the use of an appropriate parameterization, i.e. a solution of the radial equation, one can integrate the 
spin precession equations, Eq. (Tl3l) . With the solution of Eqs. pT|) - (fl6l) and ([22]) the parameter dependence of 
UN, v±, T and the tssin<I>Ar and tgcos$Ar combinations can be given. 

The second step is to substitute these quantities into Eqs. (f24f and ([25| to get the components of the N, p, 
q orthonormal triad. With the use of the transformation Eq. ([2|) one can evaluate the components of the spin 
vectors in the comoving system. 

The last step is to substitute these results into Eqs. (f27|) - (|29p and collect all the terms which belong to the 
same PN order. To get the explicit form of the polarization states quadratic or higher order spin contributions 
are neglected. The detectable wave signal arises with the use of Eq. ([1]) and the explicit form of the beam-pattern 
functions F + and F x l34l. 



4 The circular orbit case 



To investigate the effects of the orbital eccentricity on the detectable waveform we present the results of the 
method outlined above in the circular orbit limit. Furthermore, as the eccentricity decreases under radiation 
reaction in most cases the motion of the binary becomes circular at later stages of inspiral when the highest 
amplitude signals are emitted. 

To define circular orbits the decomposition §E§ of the relative velocity vector is reformulated as 
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with r = Co = 0. The connection between the above form of v and the angular variables introduced earlier is 
given as 



uo = $ + * = T 



(31) 



see Eq. With the use of these quantities the description of motion and the expressions of the N, p, q, Si 
and S 2 vectors become simpler and they can be integrated explicitly in time t. To keep the expressions of the 
polarization states as simple as possible we use the following form of v 2 given in [3U] up to 1.5 PN order 
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In this way the contributions to h + and h x can be reformulated as 
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(33) 



(34) 



The main steps of the method outlined above remain the same, only the expressions become simpler. The 
only relevant difference is that the equations of motion in this special case can be integrated in time and the 
explicit time dependence of the detectable waveform can be given. Although different post-Newtonian and 
spin-orbit contributions arise the orbital frequency, Eq. pip , remains constant (up to 2.5 PN order). As a 
consequence higher harmonics arise in the waveform compared to the Newtonian limit of this case, where the 
frequency of the emitted gravitational wave is twice the orbital frequency. 



5 Conclusions and outlook 

In this work we have discussed the orbital motion and the emitted gravitational radiation of compact binaries 
with the inclusion of the spin-orbit interaction and the relativistic post-Newtonian corrections up to 1.5 PN order. 
We have given the steps of a parameterization independent method for the calculation of the polarization states 
of the detectable gravitational waveform. All the limiting cases, i.e. binaries with one spinning component, the 
extreme mass ratio limit and the circular orbit case can be calculated from this general description. 

As a new result, the spin-orbit contributions to the equations of angular variables, including the spherical 
angles of spin vectors, are given. The usual 1PN corrections to the equations of motion are in agreement with 
previous calculations [37] . To evaluate the lowest order contributions to the waveform it is usual to substitute 
the integrated equations of motion to the quadrupole formula 34. 38J. Here we have considered the contributions 
from higher multipole moments up to 1.5 PN order with the inclusion of spin-orbit terms. 

In the post-Newtonian approximation physical systems can be studied in two steps. Since the detectable 
gravitational waveform can be unambiguously expressed as a linear combination of h+ and h x the first step 
is the construction of a method for the calculation of the polarization states which is independent from the 
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parameterization of the orbit. To do this the full description of motion, including spin precession and the 
evolution of the angular variables, is required. 

The second step determines our future plans in this area. We plan to discuss the explicit parameter depen- 
dence and the general structure of the waveform with the help of the generalized true anomaly parameterization 
[36] . To analyze the effects of the eccentricity of the orbit we will investigate the circular orbit limit, where 
the results can be expressed explicitly in time. Furthermore, with the investigation of the one spinning and 
non-rotating limits we show the importance and the details of the spin-orbit interaction. 

In the general description the next step is the evaluation of the spin-spin contributions at 2PN order, 
which provides corrections to both the dynamics of the binary and the detectable gravitational waves. This 
interaction is relevant for black hole-black hole binaries, but if one or two objects are neutron stars another 
effect, the quadrupole-monopole interaction, becomes important at 2PN order. The method presented above 
can be proven useful in the investigation of these higher order contributions. 
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